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ANALYSIS 07 SPINNING IN A MONOPLANE ffINO 
BY THS INDUCTION MBTHOD AS C0MPAH1SD WITH TEB STBIP UBTHOD* 

By L. Pqggi 



INTRODUCTION - 



It la an established fact that the phenomenon of auto- 
rotation, which causes the greater percentage of airplane 
crashes can be verified when, the wing being wholly or In 
part beyond the stall, the rolling and yawing moments set 
up as a result of a rotation about the axis of roll are such 
as to foster this rotation which does not occur when the 
wing is below the stall. 

This manifests the Importance which attaches to the 
study of forces and moments produced by the rotation of the 
airplane about an axle located In its plane of symmetry. 

The problem has been attacked by various authors, in 
particular by Fuchs and Schmidt who applied the so-called 
Btrlp method which is based upon the assumption that the 
forces and moments per unit length acting in each section 
of the wing are equal to those on an infinite cylindrical 
wing of equal section in an air flow of Intensity and direc- 
tion resulting from the apparent relative motion- of this 
section with respect to the surrounding air (reference l) . 
In other words they disregarded the induced velocities, 
which, however, are of such importance that their omission 
Is bound to result in appreciable errors, as we attempt to 
prove In this report. 



*"0alcolo dell ' autorotazlone col metodo dell 1 lnduslone in 
ala monoplana e confronto col metodo della strlscia." From 
L 1 Aerotecnica, October 1933, pp. 1255-1293. (This paper 
wbb awarded a price of 3000 lire by the Italian Research 
Council . ) 



/ 
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PART I. THEOHBTI CAL AH ALTS IS 
Hotatlon 



The notation Is prefaced by the following explanations: 

a) Whenever we refer to a quantity in which the induced 
incidence had been omitted, the symbol for the quantity la ■ 
over-scored . 

b) The subscript 0 denotes the value of the pertinent 
quantity with respect to the wing center. 

c) The signs x 1 , y 1 , K? , £, used In connection with 
vector quantities, represent the components along axis x 
parallel with the sero lift curve of the wing (constant air- 
foil section and constant angle of attack), of axis y at 
right angles to x and to the wing axis, of bxIb k, par- 
allel to the axis of motion and of axis g perpendicular 

to It and to the wing axis (fig/ 1) . 

L = wing span, 

s = ■?■—£• distance of a section of the wing from the 

3 

center (fig. l) , 

1 L = wing chord of any section: (l = nondlmensional 
coefficient) , 

7 = resultant velocity, 

ft = angular velocity of rotation of the wing, 

w = ( r— ^, ratio of velocity at wing tip with rota- 
Vox J * 

tion Cl to the component 7 0Z of the apparent 
velocity at the center of the wing. 

7 t = induced velocity. 

General Observations and Simplifications 

/71'y/r The study of an airfoil In translatory motion and at 
Jrew angle of attack which is bound up with the Bo-called 
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sBcotid-^rrorbl-em of a wing of finite span (see Pistolesi, 
Aerodinamica, pp. 250, etc.) is more complicated because: 

a) She apparent Telocity is not constant at every 
point of the wing: 

b) The angle-of-attack range in question is so large 
as to make it absolutely inadmissible, eren in first ap- 
proximation, to consider the lift curve as a straight line' 
or to substitute for each angle its tangent or sine: 

Ic) The vortices she'd by the wing are helicoldal about 
the axis of motion rather than rectilinear. 

Aside from these facts, there is the question of va- 
lidity of the well-known relation which tlos the strongth 
of the free vortices shod by tho wing to the lift distri- 
bution. And the reason for this is that this relationship 
is generally dorivod from the Eutta Joukowskl and Stoke 1 s 
theory on the conservation of the vortex theories the ap- 
plicability of which does not appear to be admissible with- 
out the other. 

On the other hand, when we consider the physical phe- 
nomenon of the formation of the vortices and bear in mind 
that it reflects the di sequilibrati on caused by the pres- 
sure differences on the top and bottom surfaces of the wing 
in accord with the change of lift across the span, it is 
seen that Prandtl's method retains its validity even for 
the case of stalling, at least in approximation, as is con- 
ceded in the present report*. 

A check of the above should be of interest and could 
be readily made in the wing tunnel by measuring the lift 
and drag on elliptical wings of varying fineness ratio at 
stalling. If this could be verified for such wings, the 
difference in fineness ratio should correspond to a differ- 
ence (constant across the span) of the apparent angle of 
attack, and would be easily computed. 



♦Various English experiment s (references 2 and 3) demon- 
strate the validity of the Eutta Joukowski theory at stall- 
ing provided that, to the extent to which .it is computed, 
the circulation cancels normally to the wake. 
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In viev of the foregoing, together with the lack of 
more reliable criteria, the proposed problem is treated by 
the conventional method deduced from Prandtl'e theory/ 

' With this* -assumption' ,it la noted that in consequence 
of consideration b) the lift and drag coefficients are no 
longer referred to the angle of attack a, a transcenden- 
tal term of the velocity component , but rather, to its tan- 

gent denoted with A. 

7 X • 

In addition, it is noted that In the expression for 
the circulation - instead of comparing the resultant veloc- 
ity 7, an irrational function of the component of motion- 
ire compare its component on axis x, which, at least, when 
disregarding the induced velocity, can, as component of ax- 
is y, have V Q ^ and co x uig as linear function of quan- 
tity" ;V ox , T oy and » x : o) y or 7 oX . 



then we write 



T = Y x I L P (x> 



"P" being a function of " A; 

which, compared with the current expression gives alone: 

00 



P = 



cos a 



The-unit wing lift and the two components-- on axis x 
and y are: 

| 7 V x I L P 
- | 7 y 7 X I L P ' (2) 



\ T x 7 X I L P 



2 



Similarly the drag per unit of wing span and' its two 
components on axi b x and y are: 
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ITT. I L S 
2 x 

\ V x V x IIS (3) 
1 V y V x I L B 

Hence 

B = Cr °° 
cos a 

This obviates the extension In the case in which the 
components of this said reaction along any other two axes 
are considered. It is noted that, so long as the induced 
velocity is disregarded, it is: 

*x = 7 0X (1 - uyft) ; 7 y = Y ox (I 0 + t» x l) 

- *o + to x ft 
= 1-^1 

Then the induced velocity is taken into account, we 
either chauge V or A. The possibility of simplifying 
the equations makeB it possible to write the circulation as 

jr = | Ct x p + r t P t »3 i (5) 

by comparing the values of 7 and A relative to the lift. 

is a function of A and of the direction of the in- 
duced velocity v A . 

' Lastly, to insure a real advantage, the coefficients 
F and B should be expressed in simple manner by means of 
the variable A. 

Then we have 

P = A x A + A s A 3 + A 3 A 3 + H a A a + H 3 A 3 + .... (6) 
B = B 0 + B x A + B 8 A 3 + . . . . (7) 
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Obviously the value of the above coefficients must he 
computed step by step on the basis of the F and H curves. 
(As to their order of magnitude, the reader is referred to 
the numerical example.) 

Regarding item c) , it is admitted that the free vorti- 
ces coil around the axis of motion, as is rigorously ob- 
served when coincident with the axis of rotation. The re- 
sults in appendix 1 then permit the substitution of recti- 
linear for helical vortices starting from the same point of 
the wins and tangent to the corresponding helical vortices. 

Determination of the Induced Telocity of the 
System of Tree Vortices 



To begin with, l.t is admitted that, in spite of the 
high angles of attack' the usual method of calculation can 
be applied to the induced angles of attack, and specifi- 
cally that, at each point of the wing, they release a vor- 
tex of intensity equal to the derivative in that point of 
the adhoring vorticity and that 1b of the circulation. 

Then, with v^ K and v^f as the induced velocity 
components along axis of motion K and along axis £ at 
right angles to it and of the axis of the wing, the approx- 
imation gives 

n 

J. 

d —j- 

d v 1K = + -i- y£ Bln(tan" X )((0 K V) (8) 

4 TT F - f I 



1 ■— i- 

d v i£ = - * — Tg t cos(tan- l )(0) K £») (9) 

1 » 4 TT I - \ 1 

i = reduced abscissa at the point in which the induced 
drag can be determined, and £' = that at the point where 
the vortex is released.* 



* "So t e : 
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low we can write* 

cob X (tan" 1 ) X x = 1 - 0.1 x - 0.02 X xa 
Hin X (tan" 1 )X x = 1.05 X.x - 0.35 X x» 
which, expressed, by a multiplicative coefficient, become 



(10) 



a r 



+ c) = 



= -drf a6 . + (ae + b) + aea + b6 + c 1 



(11) 



By integrating (ll) relative to i 1 we obtain - when 
considering that the circulation about the wing tips is 
certainly sero, otherwise we should have concentric vorti- 
ces at these points -: 

d r 



4 TT 



(1.05 i u> K - 0.35 e a u K S V! 1 fTT. YF * V " 

r 



- 0.35 & K a f* -i- L 
■ 2 



^ e 



4 TT 



d r 



(1 - 0.1 1 (0 K - 0.20 f W S K ) Z 1 ~ 

r 

-r d£ 

- 0.20 <o K a y^-jL- l 



dE' - 



♦Note the comparison between the values of (10) and the 
effective values of the particular functions in the appen- 
dix. 
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or for (10) 



4 tt 



d r 



Bin pan-^K \ ) j-i-j d £• 



d V 



~ 0.35 (V^-L L 
2 



a £ 



V i£ " 4 TT 



cob (tan T (w K 



dF 

l 1 2 1 



e - v *v 



d e« - 



- o.?o <o K a / 1 i L 

2 



d e- 



(12) 



and which is, with exception of the terms containing 

T d £' which usually are much smaller, precisely the 

resultant induced velocities obtained when the vortex fil- 
aments aro parallel to that leaving st tho particular point. 
Consequently they are perpendicular to the direction of the 
velocity in that point and their magnitude will he that of 
a system of copl&nar vortices of Identical strength. 



■analysis of Circulation Di atrl hut ion 



IFe begin with the termB containing j^j, T d £' . By 
substituting tho vector Bum 7 Q 1 for the apparent velocity 

of the center of pressure Y Q the circulation distribution 

can be computed when taking 7 0 1 as the apparent velocity 

and considering for the induced velocity only the first 
terms occurring in the second term of (12), which corre- 
sponds to. the induced velocity when the vortices are co- 
plenar to that shed at the particular point. 
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On the oasis of the defined dl stribution we can deter- 



The objectionable feature of this procedure lies in 
the difficulty of following the calculation based upon a 
sufficiently established apparent velocity, which has no 
other interest when we wish to pursue the calculations of 
the entire angle of attack range. 

.1 

The term jL x T d £ 1 is mor tractable with the strip 
method, bo that V 0 ! can be determined when the velocity 
7 0 ' 1b known. This method' involves no appreciable error, 
as the terms are generally quite small in comparison with 
the remainder of (12) . 

Putting 7 Q ! instead of V q we can determine the in- 
duced velocities on the basis of the abstracted terms with 

J_ 1 T d The following- figures always refer to a ve- 

locity at the center of the wing equal to V 0 1 . However, 
for the sake of brevity we write V Q . Then even the 

smallness of the terms enables us to differentiate between 
the two quantities, an.i we retain the notation v^ to in- 
dicate tho part of the nertinent Induced velocity. 

As to the statement about the r ectangulari ty of the 
induced velocity to the apparent velocity the induced 

v i 

angle of attack results in -y and the effective velocity 

will practically be equal in magnitude to the apparent ve- 
locity, hence the increment due to them in the circulation 




will be |v 



d 0. 



'p-oa 



d a 



Compared with (5) it is 




(13) 



which serves to calculate function P ^ . 

How we write into (5) the ? slue of the induced veloc- 
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ity derived from the preceding considerations, bo that: 
t - | C^ox (1 - w y i) PI + V 1 t 4 ] = 



"= i- V 



i , 



n (14) 



(i - o>y e) >i - 



4 IT 



d £• 



The result 1b an Integral equation which can be treat- 
ed w.ith the method' .perfectly analogous to that employed In 
the resolution of the "second problem on wings of finite 
span" cited above. 

T.or example, when GHauert'e method (described In de- 
tail In Pistolesl's treatise) which Is particularly ap- 
propriate, Is used, we can put: 



T =' 7 0X L Sa a sin n 0 



(15) 



whereby: 

cos 0= i 
Then we have (Bee above): 

. r 



(16) 



4 TT ^ 1 



2 L d p m ~ 7 ox . s n a n sin n 0 = vi (17) 



3 sin 0 i 



and (14) becomes 

_ CO 



V 0 xL 



2 a u Bin n 0 = 



= | [(1 - o) y cos 0) P I - S n a,, sin n ej 
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08 - r . Pi 1 l n i 

Z a_ Bin n 6 1 + — i— - * ? M 1 - u T co> 8) (18) 

i ^ I 4 Bin 9j 2 * 

Note: Incidentally it will be Been that this integral 
equation Is substantially the same and gives therefore the 
' same values for the circulation as that relative to the 
"second problem". for a wing of finite span moving with uni- 
form forward speed 7 0 and for which the lift coefficient 
will be proportional to the total effective angle of attack 
(apparent *Y 0 and Induced y ^) according to the ratio of 

proportionality A. The chord c and the apparent angle 
of attack 7 Q are defined from 

2. = L P 4 ' I 
Zi A 

7 Q = £-,(1 - u) y cos G) 
° P i * 

Thus it is apparent that the variation of the deriva- 
tive of the lift curve corresponds to a congruent change in 
apparent angle of attack and wing chord. 

The induced velocity then will be: 

00 

v, = - Y n * En a. sin n 8 (19) 
1 2 Bin 9 i 



Determination of Coefficients e^ 

Equation (18) being satisfied for all points of the 
axle of the wing exactly defines the infinite coefficients 
Bq. Their determination requires In general the resolution 
of a system of Infinite equations obtained opportunely by 
checking (18) for an infinite number of points on the wing 
axi s . 

In practice, since the coefficiente - being relat- 

ed to a development in Fourier series - approach zero when 
n = Infinite, and usually are small enough to be negligible, 
when starting from n = 6 to 8, we can confine ourselves to 
the first terms of the development. 
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If they are m In number the relative m coefficients 
can "bo obtained: In two ways: 

The first, which has the advantage of not requiring an 
analytical term of the curve of coefficient F and P^ , 

simply consists In prescribing that (18) be complied with 
a.t m points of the wing axis (obviously the tips must also 
bp identically satisfied) well chosen. 

» 

The second method, which we pref en, has the great ad- 
vantage of giving T and f.. (circulation with or withou't 
considering the induction effect) in the same analytical . 
form, and .which, lending Itself readily for comparison of 
both cases , is a natural extension of Hunk's method for the 
elliptical wing at low angle of. attack. 

It requires the development of 

L (1 - cos 6 ) P I 

4 sin 6 
In Fourier series, to wit: 

1 — ' ■ 00 _ 

*■ (1 - cos 6) P I = Z a sin n 9 (20) 

2 i n 



P. » I 

~ = ''n 



1" • 



= £ t_ cos n 9 (21) 



with 



4 Bin 6 1 

The first term of (18) can,, as we know, be expressed 

CD 

Z k n sin n 6 



wherein k a is a linear function of a^, 

'■■ ' To complywith (lS)it Is! - ■ 

kn = «n 
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" '\ . She result la a system of infinite equations of infi- 
nite Unknown values a^* - Then, If -the development contains 
only a limited number of terms (say, two or three), the 
problem Is markedly simplified, with the result that a num- 
ber of coefficients cancel 

Vote: In the above treatment it was tacitly assumed 
that the Increment of the lift coefficient was in every 
point proportional to the induced angle of attack. JTow, 
since the induced angles of attack can assume, as we know, ap- 
preciable values, this simplification consisting In final 
analysis, in substituting at each point of the lift curve 
the tangent of this curve, can involve considerable errors. 
Although we were not concerned with this as we did not wish 
to complicate our problem too much, it nevertheless should 
be Interesting, even for the effects on other problems, to 
see how to obtain a closer approximation. 

We have: V = velocity, a Q = apparent angle of attack 

T=Ll V C p = V L Z la sin n 6 

the value of the circulation in a certain approximation (for 
example by strip method) ; 

7l = 2 ell 9 ? ^ a n sin n 9 

the corresponding value of the induced velocity by the same 
approxlmati on ; 

T = LV Z a^ sin n 8 
the effective value of the ^circulation; 

_ V 00 

•v = 2 n a n sin n 9 

1 2 sin 8 1 

the effective value of the induced velocity. 
.We write 

£=T S„ .inn l.kttA ft- T 1> ) M 
1 1 I * do TJ 



( 
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This equation t^ similar to (14), but the . incremental . terms 
and the quantity are much reduced, hence, the approximation 
is hotter. 

Induced Velocity* 

Equation (19) yields: 

v iK = v i sin(tan *) 

1 - Hi g cos 0 

s tx. -i\ cos 9 

7j. = v* cosUan J- - 

H 1 1 - a>£ cos 9 

which may be .written in the form of: 

v ox » 

T iK = _ — fi 2t n sin * 9 
2 sin o i 

^ox oo 
S 2 sin 6 i 

The b n and c n values are easily computed, when 

boaring in mind (12) as functions of the coefficients 
(see Fart II) . 

Similar terms are obtained for the induced velocities 
projected on axiB x and y. 

Torces and Moments on the Ting 
a) Torces and Moments Due to lift. 



Tirst we abstract the induced velocities. Baaed upon 
(15) the circulation is then expressed as 



r = * ox L S an sin n 9 



•Contrary to the general rule the induced velocities are fig 
ured positive according to the positive direction of the cor 
responding axis. 
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The force components along axis K and £ due to the 
lift "per unit length of the wing are: 

* _ _ 

P K = - p "/ L V* r i s = - P T 0I : / (% oob 8 T sin6 d 9 
a 



-= - p 



TT oo 

Z 
1 



7 ° xa La o) K / Q S Sq sin n 0 cos 6 sin 6 d 6 



. = p/! ? K rds = P Vqx \ — (o £ cob 9 TsinGde (22) 

b <o O |COS Oq J 



= p L*S£ 



•^-r - /" S Iq Bin n 9 sin' 9 d9 - 



COB (X 0 0 i 



TT 



- u>£ S an 



S an sin n 9 cob 9 Bin 9 d9 



and the moment components are: 



¥ - - t„t3L 3 f^K n w 1 

M PS = " P 5 7 £ dz = " P 4 7/ 0 [«n n9sin29cos9d9| 

s a * L — 1 

I 3 3 r 

V = p £* ^Fsds = PT °^ aL3 |j o ~^/o | I n ain n 9co B 9sin9d9 - (23) 

- Z aV. Bin n 9 sin 3 9 cos 9 d 9 

3 o i ^ J 



and is 



n JO for n * 1 

/. sin n 9 sin 9d 9 =<_ 
0 ]J for n = 1 



TT 



/ sin m 9 sin n 9 cos 6 4 8' 

o 



0 for |n - m[ $ 1 



> (34) 



TT 



^ for Jn - m| = 1 
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ire finally arrive at 



! PST o 



| P S 7 0 3 



3 Tt \ W K 
= -COS a 0 HTg 

— 2 ■ 2 2 3 



= cob <io 



tt X f a x 
2 (cos a 0 



► (35) 



3-37^ = - coe 3 a 0 X 2~ [a x + a 3 ] 



J" p V o S 



I p 7/ SV 



= 00 8* a 0 



4^4 

4 ^ 



eoa.oco 



[a 



-when taking into account that L can be replaced by 
(S = area, \ = aspect ratio) and introducing the force and 
moment coefficients. 

The Induction due to (22) can be allowed for by sub- 
stituting ag for &J, , likewise V QX co K cos 6 and 



ox 



( — 1 — ojf- cos G J can bo replacod by tho moro com- 
\cos a 0 s» / 



plete equations 



_ / E c. sin n 6 

t 03C C«k co s 9 + * » 

2 sin 6 



) 



ox ( cos 6 + - — i — * S o n sin n e") 

UJS - \ cos (in t 2 sin c 1 n / 



(26) 



cob a 0 
Then we have 



f pK = - cob 



fp£ = cos' 



"o + i" f *n on] 

2 [cob b 0 2 ** 1 J 



(27) 
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and again 



.a/ r Ob -i 

mpg = - cos»<x 0 n-^ » j| WK (a 1 + a 3 ) + £ Z(an "c n+1 + c n an +x )J 



m dK = co&a 0 T— L w t (a, + a 3 ) + 

4 jjcoB a 0 2 «• 1 



1 00 "1 
+ g Z(an Cn+1 + e n a n+1 )J 



(28) 



Analogously, the comparison of (27) with (26) gives the in- 
duction effect: 

f P i ~ ? P £ = " cosa a ° [l^ (&s " S a) + | | an en] 

for the momenta. 

It will be observed that the terms containing the dif- 
ference aj - ai denote the effect of the change In circu- 
lation, whereas the termB containing b n and c n represent 

the effect - for equal circulation - of the directional 
change of the velocity. 

b) Forces and Moments Due to Drag 

On the basis of (3) we have: 

d * - \ P * Vox (l -«y*+ =~) *f ■ de 

v t n IL 8 



T ox 



Putting 



N v ox 



2 d n sin n 6 
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we similarly derive: 

1 

f 



'■ tt x r & i i°° "I 

? lrk = °o« a CL 0 — 3 — - - 0>£d a + - S d n bj 

TT \ fl 1 00 ~| 

rS = coaa a 0 — g <°k d 3 + 2 ? dn Cn 

* [cos a 0 2 

+ I f^n ^n+x + dn+x b n )J 



f 



m r £ = coe a a_ r X 3 | '■ w A (di + d 3 ) + 



m rK = cos- 



We omit the terms of the corresponding coefficients 
for the case of negligible induced velocities- since this is 
easily deduced from the above formulas. 

How let us note how by a similar method the force and 
moment coefficients with respect to axis x .and y can be 
directly established, which in special cases, such as of ro- 
tation about axis x, may make the calculations more tract- 
able. 

Wo havo with simple developments: 
f px = - coe* Oo ^ [l 0 a x + I u) x a a + || ^ c n J 

fpy = cos 3 Oo ^ L - i ^ a a + i I an b n 

\r 1 1- (3 i } 

m py = -cos a a 0 j- ^ fe U 0 a 3 + | w x(a 1 + + - S (apCn+i + cna n +i)J 
m px = cos 3 a 0 1 X^k- |* y (a 1+ a 3 ) + l| (■ ta » B+1 + 

Naturally b n and c n must be computed with respect to 
the new axes. Tho total force' and moment coefficients are 
obtained by adding tho moments due to lift and those duo to 
drag. 
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PAST II. HUMBBICAL APPLICATIOH 

The mathematical Interpretation of the effect of the 
induced velocities on the momenta acting on the wing con- 
sists in applying the results of the preceding analysis to 
the determination of the rolling and yawing moments of an 
elliptical wing of aspeot ratio 6. Shis determination ex- 
tends to various values of mean angle of attack and speed, 
of rotation for both Included and disregarded induced ve- 
locity. 

The results of these calculations being given else- 
where, we give here only a brief summary of the analytical 
method employed in the determination of the various perti- 
nent quantities. 

Characteristic- Wing Curves 



The wing characteristics are defined on the basis of 
the curve of the lift P and drag H coefficients in 
terms of variable A = tan a by moans of the following 
equations (empirically stable in following tho courso of 
the corresponding experimental v&luos relative to the wing 
taken into consideration in the work of ffuchs and Schmidt*): 

P = 5.75 A - 0.55 A 3 - 17.60 A 3 + 13.90 A* (l») 

H = 0.05 + 1.1 A* (3i) 

These curves are illustrated in figure 2, while figure 
3 gives the corresponding curves for the coefficients Cp 
and C r (with conventional notation) in comparison with the 
above experimental figures. 

As previously pointed out, the intersection of P and 
E with Cp and C r is directly given from 

_L iL i 

Cp C r cos a 
a =. absolute angle of attack 



*Ab a matter of fact, this work treats the values of Cp and 
C r for a finite aspeot ratio. However, we are not concerned 
with this except as to the trend of the curve which at least 
at high angles of attack slightly varies with the aspect ra- 
tio. 
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Tigure 3 gives onl7 the curve of the Fj/ values, 
readily obtainable from . 

Determination of a^ 

According to the •. explanation, the coefficient an are 
factors of the development of ...... 

.J- (1 - w y cos 8) P I . 
In ffourior Berles; 
or, for the elliptical wing, of 

-V (1 - u> v cos0)'Fsin8 

TT K * 

This development is readily effected when o> y = 0. In 

the opposito case it is more ezpodient to express tho vari- 
able 



. - A 0 +o) x cos 9 

1 - U)y COS 9 



in tho form of y iu x 



W y 1 ■- lD y COS 0' 

Then we substitute (1*) for F which, written in (3), 

gives 

S ^ sin n8 = [A Q + u> x cos9] U 0 + ; = j^ol g + 

+ Tl - » y cos 9) B + (1 -a) y cos 0)*J Bln 9 

where the coefficients^ A, A a -, Ag , are simply functions 
of the mean angle of attack (or what is the same, of Its 



tangent A) and' of ratio j~ 

- ■ . 7 • ' 

Thus the problem narrows' down to the development of 
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Bin 9 

i - (i) y sob e 
s£a-fi 

(1 - «y 008 6) a 

Bin e 

(1 - d)y COB 6) 3 

In Fourier Beries, which 1b easily effected. 

Pl'l 
4 Bin 6 



Approximate Evaluation of 



In the particular case of the elliptical wing, this 
becomes 

Pi' 

TT \ 

The exact evaluation of this term Involves consider- 
able complications In the calculations, at least, when we 
wish to apply the method of calculation Indicated In Part 
I, and which Is followed in the present example. 

As otherwise the torm itself is relatively small with 
respect to unity, it is added in tho first term of (18) so 
as to give it a sufficiently approximate evaluation. This 
may be done as in this particular case by constructing 
graphically by point s the curve giving this term as func- 
tion of tha. vjariable 0, 

This accomplished, the curve Is assimilated to a par- 
abola after which the said term Is readily put in the follow- 
ing form: 

Pi' 

-t- = t Q + t COB 0 + t- COB S 0 
7T \ ° 1 8 

Calculation of a- 



With due regard to the preceding approximate term, we 
can write 
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f a n ainne[l + r ^ r P.] = .... 
- ^ (1 + t 0 ) - 1 a, t a +a 3 t 1+ 2 a 3 t.J.i 

+ £a a (1+2 t 0 ) + | tl a 3 + 2 a 4 t a + 1 a x t x J 
With (16) in mind, (18) gives: 



n e + 
Bin 20 + . . . . 



C 1 + *o - ^ ) a i + *i a a + | *a»s = »i 



1 3 

5 t 1&1 + (1 + 2t 0 )a 3 + - t x a 3 + 2 t a a 4 = a a (3«) 




|- tga^ t x a a + (1 + 3 t Q ) a 3 + 2 t x a 4 + 5 tg a B = a 3 
and for the following terms in general: 

t : an-i + (1 + n t 0 ) aj + 

n + 1 a + 2 

+ — 2 — S *n+i + — 2 — *a»a+a = &n (3») 

If the consideration is limited to the first m terms 
by assuming the others to he negligible, it results in a 
system of m equations of m unknown a x a 3 a Q . 

This can be resolved, for example, hy a method of pro- 
gressive approximations the first of which may he obtained 
by assuming all terms of the principal determinant to be 
aero, excepting those of its principal diagonal which usual- 
ly have much higher coefficients than the other terms. Then 
the thus obtained valuos may bo put in the torms first as- 
sumed as zero, considering them now as shown, which yieldB 
a second approximation. The procedure can be repeated to 
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obtain various progressive approximations. In the case, 
however, where t with index greater than 1 ie eero it 
1b more expedient to use tho conventional method of resolu- 
tion by substitution. 



Calculation of Factors b n and c n 



of 



These coefficients are the results of the development 
.in (tan *)( - Jg. — a )| n ^ <lnn e 

cos (tan-^C C0B 6 - )f n a,, sin n 9 

^ 1 - cog COS 9 ' 1 

By expressing . cob — 0_ ln tne approximate form 

1 a. (I) £ CO B ^ 

<» K cob 9 +(Dk iDj cob s 6 and introducing equation (10) while 
disregarding the terms with coa 6 at powers higher than 
two, wo can write 

Bin (tan -1 ) C0B 9 = ~h, cob 6 + h a coB 8 9 

1 - u>£ cob 6 * 

cob (tan" 1 ) a * C0B 9 a = 1 + * x cob G+ kgco^e 
1 - tog cob 6 i-a 

where h and k are constants which are easily computed. 
Then we have 

• ^ = (l + fe.)^ + k iaa +|.ka a 3 



and in general: 



24 



E.A.tS.A. Technical Memo'ranaSim ITo. 747 



+ °_±J, v + n + a k 

b Q can be obtained in the same manner by substituting 
k for h and deleting the.- terms 1. .within C J (which can 
be called k Q = 1, while h 0 = 0) . 

It should be noted that the above method can be extend- 
ed to include the case where the directions in which the in- 
duced velocities are to be determined may be other than axis 
K or £, and the particular case when in direction of axis 
x and y. .... 

Coefficients d and d 

She d coefficients can be determined in the same 
manner as the a coefficients. Vow,' however, the results 
are much more simplified becasue of tho smaller number of 
terms in the £ equation compared to that for P. 

The d coefficients can be defined in the same way by 
any approximation considering only the fir at two ■ terms of 
the serieB development of the induced angles of attack so 
that the effective velocities can still be assumed as being 
linear. 

In this manner it is shown that practically any part 
included in formulas (25), (27). (28), and (29), can be cal- 
culated, which give the -forces and moments along axis K 
and £ or along x and y. It may be found expedient to 
■use the one or the other according to whether the. -calcula- 
tions insure simpler results with one than with the other. 
To illustrate: -if w y = 0, it is preferable to use the sec- 
ond, if wg = 0, to use tho first. Jrom the results of the 
first, we ban readily chango to those of tho second by a 
Bimple change of axes. 

numerical Calculation 



Case 1 .- A 0 = 0.6, rotation about axis k, which' • 
gives: 
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Ok = 1 




«x - 
u x = 



0.0833 



0.1666 



0.3600 



gftflg a- 



±0 = 0.6, 



potation about axle x(<^ = 0)* 



<°x 



= 0.10 



= 0.20 



(D 



= 0.30 



x 



Caaa 3 .- A Q = 0.3, axis of rotation and angular ve- 
locity. 

The above theory ie, strictly speaking, applicable only 
to the first case, since, in other cases, it does not satis- 
fy the system of vortices leaving the wing under the assumed 
conditions, i.e., helical vortices about the axis or rota- 
tion. However, It seemed advisable to effect the calcula- 
tions for these cases also, so as to obtain some idea of the 
degree of Influence of the induced angle of attack when the 
rotation is other than about the axis of motion. 

The following tables Bhow the results obtained in case 
1, while figures 4, 5, and 6,- give tho curves of the moment 
coefficients for the three cases. .It will be noted that, 
for case 1, these coefficients are plotted against axis K 
and g, and for cases 2 and 3, against axis x and y. 
This may seem InoonBi Btent , but it actually is not so, when 
It is observed that, in this manner, the moments about the 
axis of rotation and the rectangular axis to that of tho 
wing are considered. It affords a better characterisation 
of the wing reaction to the rotation. The curves and tables 
show for each case tho Influence of the induood angle of at- 
tack to be appreciable. It is Interesting then to note how 



♦It is easily seen that the two cases of speed or rotation, 
in the first and second cases, approximately correspond to 
the same ratio between the resultant speed of rotation and 
velocity V ox . 
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the result a of this effect In contrast to that which a pri- 
ori considerations would load us to expect. Take, for ex- 
ample, the case Ag = 0.3. Here we may be Induced to con- 
sider the lift curve as "being straight and to apply the 
methods which other authors employed from such hypotheses. 

The result would be a decreased effect in induced an- 
gle of attack with a diminished rolling moment, contrary to 
that shown in figure 5. The fact is explained by observing 
that, while in the case of absolutely straight lift curve 
the first term (constant) of the development of the induced 
velocity in Tourier series results in a lower lift coeffi- 
cient constant over the whole wing and consequently produces 
no rolling moment unless the curve is rigorously Btraight, 
but has, as in the case in question, a derivative less than 
the high angle of attack values; said decrease is less for 
the lowered wing, and produces a rolling moment of the same 
sign as that obtained by the strip method. 

Tor correct interpretation of the results hereinafter 
it Is specified that the moments opposed to rotation are 
assumed to be positive , i.e., in inverse direction of posi- 
tive rotation. 

The data in tables II and III may also appear complica- 
ted at first glance. This Is due, especially with regard 
to the rolling moments, to the fact that the induced angle 
of attack has a marked influence, and which increases as 
the moment increases when assuming the value 0.260. How- 
ever, we note that this Is due to the fact that the rolling 
moment is the algebraic sum of two quantities, lift and 
drag, which, in this case, are of opposite sign and rather 
.high value and such that, while one (lift quota) increases 
in absolute value as induction effect, the other decreases. 
In any case, the increment of the two terms is still note- 
worthy, even if taken separately. For example, the first 
increases by about 48 peroent. An explanation of this fact 
is seen in the analysis of the equation system (18) which 
"gives the coefficients in terms of a^. From this, 

bearing In mind that coefficient ta assumes high values 
(about 0.133 in this case) as can be seen from the consider- 
ations of the curve of the Pi 1 coefficients, it appears, 
for example, that a* , since it has no effect on coeffi- 
cients f and m, may have a profound effect on the coef- 
ficients a 1( a a , a 3 , and consequently on coefficients 
m and f . 
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The above Increase of the coefficient b aoeme perfectly 
natural when coordinated with other resoaroh data obtainod 
by strip method. From this (boo, for example, fig. 29, 
Luf tf ahrtf orschung, February 1929) it appears as a small 
angle of attack change which may be due to. the induced Te- 
locity raised by changes in moment coefficients of the or- 
der of those found. 

as for the rate of rotation in the pertinent field, 
there are no appreciable values at which the rolling moment 
cancels, corresponding - that is, except aa indicated - to 
the case of simple autorotation, at least, in caseB of ro- 
tation about the axis of motion or about axis x, and for 
an angle of attack in the median plane corresponding to. 
A = 0.6, we may deduce by extrapolation of the computed 
curves that this occurs at much higher rate of rotation than 
is obtained with the atrip method. 

The exact calculation of auch rates of autorotation, 
which may be UBeful for experimental verification (model in 
wind tunnel mounted free to rotate about its axis of motion) 
waa omitted here, simply because there are corresponding 
maximum values of the angle of attack greater than those at 
which the experimental curves of the wing coefficients are 
sensibly coincident with those given by the analytical ex- 
pressions used for the calculation. 

This, naturally, does not minimise the generalisation 
of the procedure, because it is always possible to give an- 
alytical exproaaiona with greater number of terms, which 
would coincide over a greator rango of angle of attack with 
tho experimental onos, or else because the analytical trans- 
lation of tho exporimontal curves la not an essential char- 
acteristic of tho method. 

However, it was not deemed necessary to effect ulterior 
calculations on the basiB of the above modifications in or- 
der to determine tho rato of aoro momont, because if we sep- 
arate the experimental verification from the above appllca*' 
tion they have not the great importance which they are sup- 
posed to have, they would not rigorously correspond to the 
true conditions of autorotation, for which, besides nullify- 
ing the rolling moments, it also requires the concurrent 
cancelling of the yawing moments. This, as may be readily 
proved, does not occur in tho conditions of flight undor 
consideration, wherein tho path of the center of gravity is 
maintained In the plane of symmetry of the airplane. The 
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rate of autorotation requires In general the consideration 
of a rate of yaw, which complicates the problem to such an 
extent as txt,.maxo_ It. almost intractable. 

TA5L1 I. 

Principal GoofflolentB for the Determination 

of the Torcea and Moments on a Wing 
Casex Rotation about Axle of Motion A = 0.6 



0) K 


0.083333 


1 

1 0.16666 


0.260 








1 w . \JO f J 


a 2 


- 0.0065 


- 0.0126 


- 0.0116 


a 3 


- 0.0001 


- 0.0005 


0.000 


a* 


0.0010 


0.0291 


0.1382 


*i 


0.0477 


0.0486 


0.0504 


d a 


0.0113 


0.0227 


0.0352 


d 3 


0.0006 


0.0015 


0.0251 


a i 


0.0784 


0.0715 


0.0772 


a a 


- 0.0084 


— 0.0486 


- 0.0178 


a 3 


- 0.0122 


- 0.0227 


- 0.0067 


a* 


+ 0.0021 


+ 0.0015 


0.0205 


*x 


0.0477 


0.0441 


0.0476 


d a 


0.0113 


0.0259 


0.0391 




0.0004 


0.0018 


0.0414 
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TU3L1 II 

Moment Coefficients Ik- and m^ 

(Induced Telocity Disregarded) 



o) K 


0.8333 


0.1666 


0.260 


a) 


Moment coefficient 
(rolling momont) 


















- U.U62U 


acq A 

- 0.0580 








+ 0.0098 


+ 0.0206 


+ 0.0356 










- 0.0414 


- 0.0224 


b) 


Moment coefficient 
(yawing moment) 


m"g 














- 0.0136 


- 0.0278 


- 0.0434 








+ 0.0559 


+ 0.1130 


+ 0.1740 








0.0433 


0.0832 


0.1386 
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TABLB III 
Moment Ooeffloients mK and mg 
(Induced Telocity Included) 



(l) K 


0.8333 


0 lfififi 


V . cow 


a) 


• 

Moment coefficient 
















- 0.0410 


- 0.0715 


- 0.0859 








0.0086 


0.0176 


0.0305 










- 0.0539 


- 0.0554 




Moment coefficient 


mg " 














- 0.0183 


- 0.0354 


- 0.0505 








0.0656 


0.1276 


0.1932 










- 0.0922 


- 0.1427 
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APPENDIX 

I. Calculation of Induced Velocities of a Vortex Tllament 
Colled Spirally Around an Axle 

We refer to the notation given In figure 7 and con- 
fine ourselves to the determination of the velocities V t " 
V p » Induced at points of axis Z parallel to axis T 
and Z. 

Consider a vortex filament of strength T which 
leaves at a point on abscissa r of the same axis. Let 
V 0 he Its rate of advance In direction of axis Z and Cl 
Its rate of rotation about the same axis. Then the vortex 

Vn 

filament has a pitch of p = 2 tt ^ and the Inclination 

of Its tangent on Z Is Gs. . 

^o 

Then we have at a point of abscissa r' ' (Plstolesl, 
Aerodlnamica. p. 472) : 



v a « = — r f° Q 2T-^a[r - r' cos 8]d 9 



4 TT 



(1) 




with 

a 

H = r» + r'a - 2 r r« cos 9 + -r§- 9 a 




By putting H In the form of 



we can write 
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in oloee approximation*. 

This expression written in (l) affords: 



v i - - — r 

a 4 TT 



CO 



(r - r») d 6 



[(, - + (£)%•] 



737 + 

.I /a 



+ r 



(1 - 00s 9 ) 



d 9 - 



- 3 r r' J 



V 



4 TT 0 



00 



00 

(l - cos 8 ) (r - r 1 co s 9 ) 

[(,-,. >• + ( £)V] 

(r» - r) d 9 

[(r - r«)« + e a ] 3/s 



d 9 



+ r/ 



or- 



1 - cos 9-9 sin 9 



, ,°°(l - cos 9 [r» - r (cos 9 + 0 sin 9fl a _ 
™ 3 r r" J _ -r ' - « d 8 

0 



On the other hand, we know that: 



r 1 - r 



d 9 = 



To r' - r 



♦The convergence of the above indicated development is in- 
sured from Doing consistently 1 - cos < 9 a so that in 

2 

the most unfavorable hypothesis tho ratio of Incremental to 



basic term ro suite in = r r 



a ratio which, as 



pointod out olsowhoro in tho roport, is conslderod < 0.25. 
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therefore we can write synthetically : 



▼ 1 

*a 



r " r - 1 v o , ~ ■ ,\i 
^ = h [m + £ u -* + * "■>] 

These expressions are useful for comparing the helical 
vortex having the same initial velocity with those of a 
straight vortex leaving at the same point of the vortex. 



These are clearly expressed in the formula 

- i r ... ,.„-^Qr_ - l r 

r • - r 



vj = = — sin (tan" 1 ) tii = - 1 — cos(tan-^2. 

4 tt v o r 1 - r 4 Tr Tir 



— _ 1 r , . « 1N flr 

v t = — — — cob (tan 1 ) — 

r ' — r 4 tt V ft 



H2 



How we shall show how the calculation of coefficients 
A and B may be carried out. To compute these integrals 
ranging between 0 and » it is advisable to divide them 
into two parts, in ono of which ranging near 0 and a fi- 
nito number suitably selected, we can develop the numerator 
in series of variable 6, while in the other, ranging be- 
tween said number and infinity, the term 



2 

can be developed in similar fashion. 



In this manner we can return to the calculation of the in- 
dicated Integrals. Obviously, in the above series develop- 
ments It suffices to consider any appropriate flnito number 
of terms. 

For example, the terms A a may be approximated at 

A a = ~/ a " 2 + 8 "" 128 + 5760 d 9 + /" =-7 d 9 - 

o r 3 [a + hS3] 3 / a 

(a + b 9) 'a 
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i [>£08e de + /«aJ a 9 T . 
+ JL JL_ a[7° £fla-S d e + Aia-J d e] 

after putting 

. a = (r - r') 3 ' b = (Jf)" 

and similarly for the other terms. 

Numerical Calculations 
Consider a wing of unit semi span, a given ratio of 

V 

= 2.5, which may "be considered as a minimum attainable 

practical limit in spinning, at least, within the range of 
validity of Prandtl's theory, and visualise a vortex leav- 
ing at the point of abscissa r = 1. We determine the 
values of the above coefficients for the points on abscissa 
0 and - 1 (the calculation of tho coofficlents B for 
tho point on abscissa 0 is omitted, sinco thoy aro ob- 
viously useless, being multiplied by the abscissa itBelf.) 

We have: 

r'l = A 1 = ■ B i = A 3 = B 3 = 

' 0" 0.05354 - 0.0280 . 

-1 0.03036 0.00120 - 0.00625 0.00236 

The following table shows the corresponding values of 

4 IT 4 tt 

T a ' -p— and v t ' in comparison with those of equa- 

tion (2): 

t , 4 TT — , 4 TT v , 4 TT — , 4 Tf j 

r» v a » -p — v a « — error $ v t « — v t « -y- error $ 

0 + 0.400 + 0.372 +7 - 0.972 - 0.930 + 4.5 

- 1 + 0.172 + 0.186 - 7.5 - 0.4990 - 0.465 + 6.7 
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Since it may "be assumed that tho percent ago at error 
e stahli ehed for the point r ' = -1 is a maximum and bo- 
Bldofl, since- the closer we appraoch the vorteX the error 
must of necessity decrease also, we wero constrained to de- 
duce that the substitution of . tho rootilinear for tho holi- 
coidal tortex does not produce, according to the effects 
-of this analysis, errors greater than those involved in the 
conventional approximation inherent to the nature of the 
problem itself. 



II. Check on the Degree of Approximation Obtained with 
a) cos (tan-i) x = 1 T 0.2 x - 0.2 X s 





b) sin (tan" 1 ) 


x = 


= 1.05 


x - 0.35 


X 3 








cos 


(tan * ; 


X 




Bin 


(tan- 1 ) a 


c 




effective 


for a) 


error f> 


effective 


for b) 


orror $ 


0.1 


0.99553 


0.988 




0.756 


0.09955 


0.1015 


+ 


1.97 


0.2 


0.98058 


0.972 




0.821 


0.19812 


0.1960 




0.06 


0.3 


0.95782 


0.952 




0.602 


0.28735 


0.2835 




1.34 


0.4 


0.92847 


0.928 




0.050 


0.37139 


0.3640 




1.99 


0.5 


0.89443 


0.900 


+ 


0.675 


0.44721 


0.4375 




2.15 


0.6 


0.85749 


0.868 




1.265 


0.51449 


0.5040 




2.04 


0.7 


0.81949 


0.832 




1.406 


0.57364 


0.5633 




1.80 


0.8 


0.78087 


0.792 




1.425 


0.62469 


0.6160 




1.39 


0.9 


0.74329 


0.748 




0.834 


0.66896 


0.6615 




1.11 


1.0 


0.70710 


0.700 




1.005 


0.70710 


0.7000 




0.71 
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Pig. 3. 



— — — — Lift and drag coefficients corresponding 

to P and R of Fig. 2. 
— - — — Bsperlnontal coefficients on Junkers 045 

wing. 
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1.0 
0.9 
0.8 

p 
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0.5 
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tfigure 4a.- lionents due to rotation about axLe of 
notion, co £ a 0 for Aq b 0.6 



Coefficients vt» • and n 
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ligore 4b. - Moments due to rotations about axis of 
notionui g ■ 0 for Aq = 0.6 



